Incommensurate spin-density wave and magnetic lock-in transition in CaFe 4 As 3 
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The magnetic structure for the newly discovered iron-arsenide compound CaFe4As3 has been 
studied by neutron powder diffraction. Long-range magnetic order is detected below 85K, with an 
incommensurate modulation described by the propagation vector k=(0,<5,0), 8 ~ 0.39. Below ~ 
25K, our measurements detect a first-order phase transition where 5 locks into the commensurate 
value |. A model of the magnetic structure is proposed for both temperature regimes, based on 
Rietveld refinements of the powder data and symmetry considerations. The structures correspond 
to longitudinal spin-density-waves with magnetic moments directed along the 6-axis. A Landau 
analysis captures the change in thermodynamic quantities observed at the two magnetic transitions, 
in particular the drop in resistivity at the lock-in transition. 



PACS numbers: 75.25.+Z, 75.47.Lx, 61.12.-q, 75.10.-b 

The recent discovery of superconductivity at 26K in 
iron-based layered compounds La[Oi_ x Fa;]FeAs |l| has 
triggered a new wave of investigations in the field of su- 
perconductivity. There are many similarities with the 
High Temperature Cuprate superconductors (HTC) with 
the presence of square planar superconducting layers, in 
this case derived from tetrahedrally-coordinated Fe2As2 
units, separated by charge reservoirs. All the undoped 
non-superconducting systems, including the so-called 122 
compounds AFe2As2 (where A=Ba,Sr,Ca), undergo a 
tetragonal to orthorhombic transition accompanied by 
a commensurate antiferromagnetic (AFM) order with an 
ordered moment less than 1 \xb [2j|- Chemical doping 
suppresses the structural transition and superconductiv- 
ity is observed at critical temperatures as high as 55K 
and 38K in optimally doped SmFeAs(Oi_ x F a; ) [3| and 
Bao.6Ko.4Fe2As2 Q respectively. 

Since the Fe2As2 layer is the key structural motif in 
the new iron arsenide superconductors, the ability to 
explore different topologies that include such building 
blocks is of paramount importance. The recently dis- 
covered CaFe4 AS3 [H| compound is an exciting member of 
this family, where the crystal structure is made of in- 
finite Fe2As2 layers along b but finite in a perpendicu- 
lar direction. These Fe2As2 ribbons, interconnected by 
five-coordinated iron atoms, are arranged in a rectan- 
gular cross pattern. Magnetic susceptibility measure- 
ments IH reveal the existence of two transitions at 
T/vi ~90K and Tjv2 ~26K, with antiferromagnetic cor- 
relations along the 6-axis at Tni- The heat capacity mea- 
surements show a A- type anomaly only at Tjvi, and an 
extremely weak change of slope at Tat2@- On the other 
hand, the resistivity measurements show a weak anomaly 
at Tjvi, but a more pronounced drop at Tn2 although the 



magnitude and direction differ in the two studies 0, @]. 
Even though such bulk measurements clearly point to a 
strong coupling between the magnetic order parameter 
and resistivity, to date there is no scattering studies of 
the magnetic state in this system and its dependence on 
temperature. 

In this communication, we report on the magnetic 
structure of CaFe4As3 probed by neutron diffraction ex- 
periments on polycrystalline samples. The transition at 
Tni corresponds to an incommensurate (ICM) longitudi- 
nal spin-density wave (SDW) along the 6-axis with prop- 
agation vector k=(0,<5,0) and 0.375< S < 0.39. The tran- 
sition is second-order with the magnetic order parame- 
ter transforming as one of the irreducible representation 
(IR) of the Prima space group. At Tn2, the structure 
locks into a commensurate (CM) state with £=| and 
the moments remain predominantly aligned along the b- 
axis. The CM-ICM transition is first order and the CM 
state is stable down to the lowest temperature measured 
of 1.5K. A phenomenological description of the magnetic 
phase transition from the paramagnetic state to the ICM 
and CM state is able to reproduce all the features found 
in the temperature dependence of physical quantities. In 
particular, the sudden drop of electrical resistivity at Tn2 
associated with a surprisingly weak anomaly in the heat- 
capacity are both explained by the presence of a sixteenth 
degree invariant that stabilizes the CM state. 

A 0.8g sample of crushed single crystals grown in a 
Sn flux [5| was loaded into a standard 6mm vanadium 
can and cooled inside a cryostat for neutron experi- 
ments performed at the ISIS pulsed neutron source of 
the Rutherford Appleton Laboratory, U.K.. The temper- 
ature dependence of the crystal and magnetic parameters 
was followed using the General Materials Diffractometer 
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(GEM). Data were collected on cooling between 130K 
and 5K in 5K steps. Higher-resolution data at low-Q 
(scattering vector) were collected on the WISH diffrac- 
tometer on the ISIS Target Station 2 at 1.5K, 30K and 
100K. Rietveld refinements were carried-out with the pro- 
gram FullProf Q • Symmetry (representation) analysis is 
presented using Kovalev's notation 

The data obtained at 100K in the paramagnetic 



Site i Pos. 
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Mi (/i S ) $i 


1.5K 


Fei x=0.0210(5) 

z=0.3135(4) 
Fe 2 x=0.0671(4) 

z=0.5376(5) 
Fe 3 x=0.3037(4) 

z=0. 1243(4) 
Fe 4 x=0.3187(4) 

z=0.7233(4) 


1.4(11) 
1.61(14) 0.14(3) 
1.67(20) 0.45(3) 
1.84(10) 0.01(4) 


2.14(13) 
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1.94(10) 0.10(4) 
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TABLE I: Top) Amplitude of the magnetic modulations and 
relative phases for the four inequivalent Fe sites (i=l,4) ex- 
tracted from the Rietveld refinements in the ICM phase 
at 30K and CM phase at 1.5K (see text for details of 
the parametrization) . Bottom) Matrix representative of the 
complete irreducible representation n for generators of the 
Primal' space group. The symmetry operation are listed in 
the Seitz notation. Matrices for the time-reversal operator (- 
Identity) and lattice translations along a and c (Identity) are 
not shown. €=e~ l7TS and e* denotes the complex conjugate. 

regime confirms the crystal structure previously pub- 
lished for the CaFe4As3 (orthorhombic space group 
Pnma with a=11.5844 A, b=3.7368 A and c=11.5833 
A). There are four independent Fe sites in the crystallo- 
graphic unit-cell, all of multiplicity four and in the (x, j,z) 
positions and labelled according to Tab. fl] A small 
amount of Sn impurity from the flux is also present in 
the studied sample, and has been treated as a secondary 
phase in the Rietveld refinements. On cooling below 
Tjvi «85K as shown in Fig. Q], additional Bragg peaks 
are observed with intensities that decrease rapidly as a 
function of Q, indicative of long-range magnetic ordering. 
Based on a grid search 0] in the first Brillouin zone, these 
peaks can be indexed only with an ICM propagation vec- 
tor k= (0,5,0) with 8 close to |. The three most intense 
peaks are labeled with their reciprocal positions in the 
thermodiffractograms in Fig. [TJ The transition from the 
paramagnetic phase to this ICM phase is second-order. 
The value of 8 increases on cooling to reach a maxi- 
mum close to 0.39 just above Tn 2 «25K, and changes 
abruptly below this temperature to lock at the commen- 
surate value £>=§ within the error of our measurements 
(Fig. Q]). Subsequent measurements with smaller tem- 
perature steps close to Tjv2, not shown, indicate that 
the second transition is first order, with a coexistence of 



the ICM and CM structures over a temperature regime 
of 3K. This sequence of magnetic transitions rules out 
the possible ferromagnetic state investigated by Density 
Functional Theory Lattice strictions are also observed 
at the two magnetic transitions, as shown from the tem- 
perature dependence of b showing a negative thermal ex- 
pansion below Tjvi and a small contraction of c below 

Tat2- 

Models for the CM and ICM magnetic structures 
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FIG. 1: (Color online) a) Neutron thermodiffractograms of 
CaFe4As3 collected on the low angle bank (scattering angle 
26» = 17.98°) on the GEM diffractometer. The intensity is 
color-coded, b) temperature evolution of the component 8 of 
the magnetic propagation vector k =(0,5,0). The solid line 
is a guide to the eye c) temperature evolution of the lattice 
parameters b (closed circles) and c (open circles). Dashed 
lines indicate the transitions Tjvi (panels a and c) and Tjv2 
(all panels). 

have been derived from Rietveld refinements at 1.5K and 
30K, respectively. Due to the large number of indepen- 
dent sites and the limited number of magnetic Bragg 
peaks, fully unconstrained model even with global op- 
timization procedures do not lead to satisfactory solu- 
tions. However, symmetry analysis allows us to reduce 
greatly the number of possibilities. Since the first transi- 
tion is second-order, the magnetic order parameter must 
transform as one of the four one-dimensional irreducible 
representations of the Pnma little group with k= (0,8,0), 
labeled t\ to t±. The systematic absence of the (0,k±<5,0) 
reflections indicate that the moments are aligned along 
the &-axis. This is in agreement with the decrease of the 
magnetic susceptibility along & at Tjvi- 

Within those constraints, we found that only mag- 



3 




2 4 6 8 10 

d- spacing (A) 



FIG. 2: Rietveld refinements of CaFe 4 As 3 at 1.5 K in the CM 
phase(top), at 30K in the ICM phase(middle) and 100K in the 
paramagnetic phase(bottom). On each plot, the gray crosses 
and solid line represent the experimental data points and cal- 
culated diffraction pattern, respectively and the difference is 
shown at the bottom as a solid line. Rows of markers indicate, 
from top to bottom, the positions of the nuclear reflections 
for CaFe4As3 and Sn (impurity from flux) for all tempera- 
tures and the magnetic reflections for CaFe4As3 in the CM 
and ICM phase. 



netic modes that transform as t\ can explain the observed 
diffraction pattern. In this model, the magnetic structure 
corresponds to a longitudinal SDW, as postulated in Ref. 
Q. Each Fe site (2=1,4) in position (x,i,z) listed in Tab. 
U generates three additional Fe positions by symmetry 
operations of the group with coordinate triplets: (-x,f ,- 
z),(x+i,i,-z+i),(-x+i,|,z+i). The magnetic moment 

on each Fe site (i=l,4) in the cell Rl is written m,=M, j 
cos(27r(kRi+$ i )) where j is a unit-vector along b, Mi is 
the amplitude of the modulation and $i is the phase, in 
units of 2ir. The n symmetry constraint requires that the 
magnetic parameters (M,*,$j) transform as (M i; $j + 4)> 
(-Mj,$i), (-Mj,$j + 4) for these three positions respec- 
tively. If no additional physical constraints are imposed, 
i.e. the amplitude of the modulation is allowed to be dif- 
ferent on the inequivalent sites, the magnetic structure is 
fully described by seven parameters: the four amplitudes 
Mi and three relative phases, since one of the phase can 
be fixed to zero (equivalent to an absolute phase for the 
entire structure). A series of simulated annealing runs 
with these seven free parameters, using integrated inten- 
sities of the magnetic peaks extracted from a Lebail fit, 
systematically converged to a unique solution. The re- 
sult was introduced in the Rietveld refinement leading 
to good agreement with the data and magnetic structure 
factors of 6.1% and 10.2% at 1.5K and 30K respectively. 
The value of each refined parameter is shown in Tab. 



U and the Rietveld refinements are displayed in Fig. [2] 
We could not exclude the presence of a weaker secondary 
mode at the ICM-CM transition, since this transition is 
first-order, but we could not find evidence of additional 
components below Tjv2 within the statistics and intrinsic 
limitation of the powder data. In this respect, subsequent 
work on single crystal diffraction and neutron polarime- 
try is needed to pin down possible extra modes. 

All magnetic sites have an ordered magnetic com- 
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FIG. 3: (Color online) Magnetic configuration in the CM 
phase of CaFe4As3 at 1.5K along the b-axis. For simplic- 
ity, only one Fe2As2 ribbon is shown in blue (dark grey) with 
capping sites in orange (light grey) and the projected view is 
perpendicular to the chosen Fe2As2 layer. The + and - signs 
indicate if the Fei ion is located above or below the layer. The 
Fe sites are labelled 1 to 4 according to the atomic positions 
given in Tab. [I] The As atoms sublattice is not shown for clar- 
ity. Straight solid lines mark the connectivity in the Fe2As2 
layer. Three inequivalent magnetic exchange paths between 
five-coordinated Fei and atoms in the Fe2As2 layer are shown 
respectively as short-dashed, long-dashed and curled lines. 

ponent in the ICM phase with moments smaller than 
2/.ib- The most remarkable change at the ICM-CM tran- 
sition, in addition to the lock-in of the wave- vector, is an 
increase of the amplitude on the Fei site, i.e. the five- 
coordinated site connecting adjacent Fe2As2 ribbons. It 
is however impossible to determine whether the small 
change of relative phases detected at the CM transition 
is real, for which precise refinements of the magnetic pa- 
rameters on single crystal data will be necessary. At 
1.5K, the modulation amplitude is the largest on Fei, and 
the lowest on Fe2. This is also what was found by Den- 
sity Functional Theory (DFT) calculations [f|, although 
they were performed assuming a ferromagnetic state. In- 
deed, the ordered moment on the five-coordinated Fe was 
calculated to be 2.06 hb and the smallest moment was 
lfiBOn Fe2 (labelled Fe3 in Q). Mossbauer spectroscopy 
also suggests that out of the four iron sites, only the 
five-coordinated site has the formal oxidation state Fe + , 
which is consistent with the higher moment. The mag- 
netic arrangement is shown in Fig. [3] for a single Fe2As2 
ribbon capped by first-neighbour pyramidal Fei sites. 

The origin of the modulation can be understood within 
a localized picture by considering and known magnetic 
configuration of the parent 122 compound and its analogs 
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[U, H|. In such systems, ferromagnetic rows in one direc- 
tion of the square lattice are coupled antiferromagneti- 
cally along the other direction, leading to anti-parallel 
spins along the square diagonals (AFM next-nearest ex- 
change). In CaFe4As3 additional exchange couplings are 
introduced between the layers and Fei sites. There are 
three different exchange paths marked by different sym- 
bols in Fig. [3] Importantly, the exchange integrals be- 
tween Fei and Fe2 along the 6-axis are equivalent by 
symmetry (left and right) due to the mirror plane per- 
pendicular to b. This implies that whatever the sign 
of this exchange interaction, it will compete with the 
in-plane next-nearest neighbor exchange. This is topo- 
logically equivalent to a chain with competing first- and 
second-neighbor interactions and can explain the long- 
wavelength modulation along the 6-axis. Along the other 
diagonal, competition is also present but the two ex- 
change integrals are not equivalent by symmetry. Al- 
though there is a modulation along this diagonal inside 
a single Fe2As2 ribbon (Fig. [3]), the presence of four 
ribbons per unit-cell interconnected through Fei sites on 
a non-frustrated cross-pattern lattice, does not give rise 
overall to a modulation along a or c. 

Finally, with the knowledge of the magnetic phase 
transitions, we suggest possible mechanisms for the ob- 
served anomalies in the heat capacity and resistivity by 
applying Landau theory. Taking into account the sym- 
metry of the two-component magnetic order-parameter rj 
= | r\ | e lip and rj *= \ rj \ c~ lip transforming as one 
gets the Landau free-energy density 

$(T, | V \,<p) = $o(T) + ^ 2 + ^ 4 + ^ 16 c OS (16^) (1) 

where 71 =0 at the transition to the ICM phase, which 
displays the mmml ' magnetic point-group symmetry, 
and 71 ^ at the transition to the CM phase. Min- 
imization of $ yields three possible stable structures 
for the CM phase, having the respective symme- 
tries P112i/a(ip=mr/S), P2 1 ma(ip={2n+l)ir/m) and 
Plla{tf 7^n7r/16) , with an eightfold magnetic unit-cell. 
From the refinement of the crystal structure, there 
is no indication of peak broadening at the CM-ICM 
phase, which suggests that the CM phase retains an 
orthorhombic symmetry (P2\ma). Anyhow, for all the 
aforementioned symmetries, the sixteenth degree term 
of the lock-in invariant, corresponding to k=(0,|,0), is 
key to explain the coupling to the physical properties 
and describe the unusual temperature dependences of 
the specific heat C, which shows no noticeable anomaly 
at and the resistivity p, which exhibits a sharp 

drop below Tn2- Since the 71 invariant induces in 
the CM phase a higher order perturbation A77 2 sw ji 
(Tat 2 -T)~ with respect to its equilibrium value ri 2=9 p~ 

(Tjvi-T) in the ICM phase, the specific heat C = -T§f£ 
shows, on cooling below Tjv2, a slight drop in the 



slope AC oc T(T/V2 — T)~r with respect to its linear 
temperature dependence AC = above Tjv2, as 

reported experimentally. The resistivity contribution to 
the free-energy $ p =p,r] 2 p+^p 2 gives in the ICM phase 
P e ( T )=- i uV 2 = -^(Tjvi-T) corresponding to a linear 
decrease of p e (T). In the CM phase, the coupling of p 
to the additional sixteenth-degree lock-in term yields 
P e ( T )= -ffj [A 1 (TjV2-T)+^'(TjV2-T) 8 ]. This corresponds 
to a strong drop of /? e (T) below T7V2 assuming a large 
positive value of the coupling coefficient p\ 

In summary, we have presented for the first time 
evidence for two magnetic transitions in CaFe4As3 
corresponding to longitudinal spin-density-waves, 
incommensurate below Tjvi=85K and locked-in at 
k=(0,|,0) below T A r 2 =25K, as determined by neutron 
powder diffraction experiments. Models for the mag- 
netic structures in both phases have been derived from 
symmetry considerations and Rietveld refinements of the 
neutron data. The lock-in transition can be explained 
phenomenologically by the presence of a sixteenth 
degree invariant in the thermodynamic potential that 
reproduces the anomaly observed in the temperature 
dependence of the electrical resistivity. An in-depth 
analysis of the microscopic origin of these two magnetic 
states would require first principle calculations, which 
could also indicate the conditions under which supercon- 
ductivity could be introduced in this compound. 
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